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At Auger energies only model predictions enable us to extract primary cosmic ray 
features. The simulation of the shower evolution depends sensitively on the first few 
interactions, necessarily related to the quality of our understanding of high energy 
hadronic collisions. Distortions of the standard "soft semi-hard" scenario include novel 
large compact dimensions and a string or quantum gravity scale not far above the 
electroweak scale. Naively, the additional degrees of freedom yield unification of all 
forces in the TeV range. In this article we study the infiuence of such precocious 
unification during atmospheric cascade developments by analyzing the most relevant 
observables in proton induced showers. 
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1 Introduction 



Very recently, it has become evident that a promising route towards reconciling 
the apparent mismatch of the fundamental scales of particle physics and gravity is to 
modify the short distance behavior of gravity at scales much larger than the Planck 
length. Such modification can be most simply achieved by introducing extra dimensions 
(generally thought to be curled-up) in the sub-millimiter range 0, Within this 
framework the fundamental scale of gravity M* can be lowered all the way to O (TeV) , 
and the observed Planck scale turns out to be just an effective scale valid for energies 
below the mass of Kaluza-Klein (KK) excitations . Clearly, while the gravitational 
force has not been directly measured far below the millimeter range |Q, Standard 
Model (SM) interactions have been investigated well below this scale. Therefore, if 
large extra dimensions really exist, one needs some mechanism to prevent SM particles 
from feeling those extra dimensions. Remarkably, there are several possibilities to 
confine SM fields (and even gravity) to a 4 dimensional subspace (referred to as a 
brane- world) within the (4 + n) dimensional spacetime 0. While the phenomenology 
of n large compact dimensions and TeV scale strings is very exciting on its own P], at 
the same time it opens up new scenarii in which to explore "exotic" KK-cosmologies 
0, as well as extraordinary astrophysical effects Naturally, an intense activity to 
assess its experimental validity in collider experiments is currently underway 0. 

The extremely high center-of-mass (cm.) energies attained in cosmic ray collisions 
at the top of the atmosphere are well above those necessary to excite the hypothetical 
KK modes which would refiect a change in spacetime dimensionality. Therefore, a 
natural question to ask is wether KK excitations could have a direct influence in the 
development of extensive air showers. Cascades initiated by neutrinos (with cross 
sections reaching typical hadronic values) have been extensively discussed elsewhere 



TO, 11, 12, 13, M. We concentrate here on proton-induced showers. Before going into 



the technical details of the shower simulations, let us summarize the pecularities of 
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gravity with n > that can ruffle the standard "soft semi-hard" scenario. 



Back of the Envelope Insights into KK-Modes 
Phenomenology 



In the canonical example of spacetime is a direct product of ordinary four- 
dimensional spacetime and a (flat) spatial n-torus with circumferences Lj = 2'7rrj {i = 
l,...,n), generally of common linear size = Vc- As mentioned above, SM fields 
cannot propagate freely in the extra dimensions without conflict with observations. 
This is avoided by trapping the flelds to a thin shell of thickness 6 ~ |T5[. The 



only particles propagating in the (4-|-n) dimensional bulk are the (4-|-n) gravitons. 
Because of the compactiflcation, the extra n components of the graviton momenta are 
quantized 

ki = = —, i = l,...,n. (1) 

Thus, taking into account the degeneracy on li G Z, the graviton looks like a massive 
KK state with mass 



1/2 



It is important to stress that the graviton's self-interactions must conserve both ordi- 
nary 4-momenta and KK momentum components, whereas SM flelds (that break trans- 
lational invariance) do not have well deflned KK momenta in the bulk for i/r^. < Mg 
T6| . Therefore, interactions of gravitons with SM particles do not conserve KK mo- 



mentum components. 

Applying Gauss' law at r <^ and r ^ r^, it is easily seen that the Planck scale of 
the four dimensional world is related to that of a higher dimensional space-time simply 
by a volume factor, 

*'"'V'" 1 - 2.0 X 10- f I^U^^fiV'" cm. (3) 



so that M^, can range from ~ TeV to Mpi = 10^^ GeV, for < 1 mm and n > 2. For 
n < 6, the mass splitting. 
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is so small that the sum over the tower of KK states can be replaced by a continuous 
integration. Then the number of modes between |£| and \l\ + di reads, 



dN = dlidl2 ...din = Sn^i W^dl, (5) 

where 

^"-^ = rW2) 

is the surface of a unit-radius sphere in n dimensions. Now using Eqs. @ and @, Eq. 
(^) can be re-written as 



dN = Sn-i[^\ — m"-^ dm. (7) 

From the 4-dimensional viewpoint the graviton interaction vertex is suppressed by Mpi. 
Roughly speaking, am oc M~^. Now, introducing dam/dt, the differential cross section 
for producing a single mode of mass m, one can write down the differential cross section 
for inclusive graviton production 

d^a _5„_^^^V^^n-iC^^. 



dtdm \M,JM^ dt ' 

or else, the branching ratio for emitting any one of the available gravitons 

r. - (9) 

where s^^^ is the cm. energy available for graviton-KK emission. All in all, one can 
see by inspection of Eq. that the enormous number of accessible KK-states can 
compensate for the Mpj factor in the scattering amplitude. 

Having outlined the general ideas for the production of KK gravitons in very high 
energy collisions, let us consider now the hadronic scattering of two SM particles. To 
illustrate the effect of extra dimension gravity, we will estimate the effects of exchang- 
ing a tower of KK gravitons between the hadrons, rather than the production of soft 
gravitons. As usual, the parton evolution of interacting hadrons a and b must be 
separated into: (i) the non-perturbative soft cascades, characterized by a small mo- 
mentum transfer qt < qo ^ 2 GeV and described by soft Pomeron exchange, (ii) the 
hard cascades, qt > qo, that should be described perturbatively [0. If one envisions 
a scattering process considering the exchange of gravitons, as a qualitative assessment 
the overall shape of the cross section could be written as fL^ 



^ I „4— dim /I n^ 

(^tot = cr + a [lU) 
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where cr^^ denotes the contribution from the virtual graviton exchange, and 



a 



4— dim/ 



ah 



^ab 



Here, xTh i^jb) stands for the soft eikonal defined by pO 
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^ab 



4 y?2 

^^ah 



(12) 



where h is the impact parameter, y = In s, A = ap(0) — 1, and i?^^ = R^^ + Rl + a'p{^)y. 
The parameters of the Pomeron trajectory (A and ap(0)) as well as those describing 
the Pomeron-hadron vertices (7 and R^) are set to their values in qgsjet in the 
air shower simulation [ITHi. The semi-hard interaction is treated as the soft Pomeron 



emission (soft pre-evolution) followed by the hard interaction of partons 



^^■•'^(5,6) 



I dyr I dy^ x:f(e^"^"^&) ^hard(e 



y-Va-Vb 



(13) 



where 2/1(2) are the rapidities of the Pomeron end, (Thard is the parton interaction cross 
section, is an adjustable parameter associated with parton density and Cah is the 



211]. The latter is also fixed to the value of QGSJET in 



shower enhancement coefficient 
the simulations. 

A complete theory of massive KK graviton modes is not yet available, making it 
impossible to know the exact cross section at asymptotic energies. A simple Born 
approximation to the elastic cross section leads, without modification, to a 



KK 



\L4^ . Unmodified, this behavior by itself eventually violates unitarity. This may be 
seen either by examining the partial waves of this amplitude, or by noting the high 
energy Regge behavior of an amplitude with exchange of the graviton spin-2 Regge 
pole: with intercept a{0) = 2, the elastic cross section 



da \AR{s,t)\^ 
dt s2 



~ s 



2a(0)- 



2 C.2 
~ S , 



whereas the total cross section 



^KK 



Im[A^(0)] 



^a(0)-l 



(14) 



(15) 



so that eventually > a^^. Eikonal unitarization schemes modify these behaviors: 
in the case of the tree amplitudes the resulting (unitarized) cross section 

\n\s/so) m. How- 



KK 



whereas for the single Regge pole exchange amplitude, a 
ever, the Regge picture of graviton exchange is not yet entirely established: both the 
(apparently) increasing dominance assumed by successive Regge cuts due to multiple 
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Regge pole exchange [10, T7|, as well as the presence of the zero mass graviton can in- 
troduce considerable uncertainty in the eventual energy behavior of the cross section. 
Hereafter, we work within the unitarization framework and adopt as our cross section 
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Atts 



10 



-28 



TeV 



E 



1019 eV 



cm 



(16) 



3 Air Shower Simulations 



The experimental information obtained at ground level is only indirectly connected 
to the first few generations of hadrons. Consequently, the study of the influence of KK- 
modes on hadronic interactions with cm. energies s^^^ > 100 TeV, requires correctly 
simulating the intrinsic fluctuations in the air showers. 

Let us first discuss in a very general way the possible effects introduced by virtual 
graviton exchange. The survival probability N at atmospheric depth X of a particle a 
with mean free path 

Aa = (17) 
^'"a— air 

is given by 

N{X) = e-^/^% (18) 

where mair is the mass of an average atom of air |2^, and the cross sections o"a-air 
inferred from Eq. ( |16|) are shown in Fig. 1. It is straightforward to see that the total 
thickness of the atmosphere corresponds to more than 20 hadronic interaction lengths, 
depending on the primary zenith angle. The key feature in the evolution of the shower 
is the branching between decay and interaction of secondary hadrons along their path 
in the atmosphere. The latter strongly depends both on particle energy and target 
density. 

Because of the low air density at the top of the atmosphere the point of the first 
interaction fluctuates considerably from shower to shower. However, KK-graviton ex- 
change significantly reduces the nucleon attenuation length, e.g., at 3x 10^° eV, A^ ~ 41 
g/cm^, whereas X^'^"'^ ~ 38 g/cm^. Moreover, tiny deviations on the mean free path of 
non-leading secondaries yield a small change in the shower interaction length. Namely, 
the survival probability of a secondary pion (say E = 5 x 10^^ gV) at X = 40 g/cm^ is 
reduced from 43% to 41%, and that of a kaon with the same energy from 30% to 29%. 
Therefore, one can - perhaps naively - state that phenomenological models consider- 
ing the virtual exchange of graviton towers would trigger, on average, earlier shower 
developments than a naked "soft semi-hard" scenario. 
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Test simulations runs of giant air shower evolution have been performed, choosing 
typical parameters for the experimental situation at the Fly's Eye and Auger experi- 
ments 1^^. The algorithms of aires (version 2.1.1) were slightly modified so as to 



track the particles in the atmosphere via the standard 8 parameter function, 

where u = \n.E [GeV] and the coefficients are listed in Table 1. The hadronization 
algorithm that translates the parton strings produced during the scattering process 
into ordinary particles, remains the same. 

In the simulation, several sets of protons with £" = 3 x 10^° eV were injected at 100 
km above sea level (a.s.l.). The sample was uniformly spread in the interval of 0° to 
50° zenith angle at the top of the atmosphere. All shower particles with energies above 
the following thresholds were tracked: 750 keV for gammas, 900 keV for electrons 
and positrons, 10 MeV for muons, 60 MeV for mesons and 120 MeV for nucleons. The 
results of these simulations were processed with the help of the aires analysis package. 

The atmospheric depth X^ax at which the shower reaches its maximum number 
of secondary particles is the standard observable to describe the speed of the shower 
development. The charged multiplicity, essentially electrons and positrons, is used to 
determine the number of charged particles and the location of the shower maximum 
by means of 4-parameter fits to the Gaisser-Hillas function p6[] . 



N^\X) = N^^ { / ^° exp ~ ^ ^ X > Xo, (20) 



where Xmax, X^^, A, and Xq are the free parameters to be adjusted p^. Shown 
in Fig. 2 are the resultant X^ax distributions of proton showers with 3 x 10^° eV 
and primary zenith angle 43.9°.[] The tails (Xmax > 900 g/cm^) of these distributions 
were fitted with exponentials (a e"'^"^""''''), fioating both the normalisation a and the 
exponent in the fit. The resulting parameters are: (3 = 2.6 ± 0.1 x 10~^ cm^/g for the 
4-dimensional case, and /3 = 2.9 ± 0.1 x 10~^ cm^/g for the (4 + n)-dimensional case. 
A statistically significant difference between the two approaches arises in the tail of 
the distribution. This is because the depth of such penetrating showers increasingly 
reflects that of the flrst interaction [^]. Results of the flts to the X^ax distributions 
generated by applying progressively less restricted data cuts (distances near the peak) 
lead to exponential slopes that within the error are consistent with one another. 

In Fig. 3 we show the longitudinal developments of proton showers superimposed 
over the experimental data of the world's highest energy cosmic ray shower observed 



^This is the primary zenith angle of the Fly's Eye event |28 
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Figure 1: Inelastic cross sections as a function of the cm. energy. The sohd hne 
stands for the usual 4-dimensional cross section of QGSJET, whereas the dashed line 
represents corrections coming from the virtual graviton exchange. We also show in the 
figure experimental points of the inelastic p-air cross section as observed by different 
cosmic ray experiments [BH . 



Table 1: Coefficients for mean free path parametrization, = 1 TeV 



particle 


Pi 


P2 


^3 






P6 


Pj 


Ps 


nucleons 


-59.852 


-1916.4 


-25.508 


3.2875 


925.76 


69.860 


-0.089103 


-0.12169 


pions 


-70.680 


-1500.1 


-26.015 


2.7753 


585.44 


69.425 


0.36761 


-0.12197 


kaons 


-84.984 


-953.40 


-23.677 


2.0865 


262.26 


54.498 


0.70365 


-0.10659 
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Figure 2: Distributions of X„ 



to date We selected from our shower sample those with a primary zenith angle of 
43.9°, setting the observation level at 850 m a.s.l and with geomagnetic field specific 
for the Fly's Eye site. Although at the same total energy a shower that takes into 
account the virtual graviton exchange develops faster than that modelled with unmod- 
ified QGSJET, as expected from our previous analysis, the differences in the position 
of Xmax fall within the errors. However, there are visible deviations in the evolution 
of the charged multiplicity. To estimate the amount of departure from the standard 
4-dimensional scenario we analyzed the data by means of a t^st [0. We assume 
that the set of measured values by Fly's Eye are uncorrelated (any depth measurement 
is independent of any other), and make use of the quantity 



X 



a,- 



at 



(21) 



where q is the total number of points in the analysis, a^j is the error on the Xjth 
coordinate, Xj is the measured value of the coordinate, and aj the (hypothetical) true 
value of the coordinate. The obtained results are X4/DOF = 324.89/12, x^^^^^/DOF = 
200.52/12. If in the future the situation should arise that one can be confident that 
the hadronic interactions are correctly modeled, then it will be necessary to carry 
out a more sophisticated statistical analysis which, for example, accounts for the non- 
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Figure 3: Atmospheric cascade development of proton showers [E = 3 x 10^° eV), 
superimposed over the Fly's Eye data. The error bars in the simulated curves indicate 
RMS fluctuations of the means. 
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Table 2: Particle densities [m The errors indicate the RMS fluctuations. 



R = 1000 m 
82.38 ±26.10 X 10° 
19.71 ±01.32 X 10° 
62.46 ± 26.00 x 10° 
91.73 ± 10.60 X 10^ 

R = 1000 m 
10.99 ±02.08 X 10^ 
21.61 ±01.63 X 10° 
88.28 ±21.10 X 10° 
10.40 ±01.42 X 10^ 



Gaussian distributions. 

It is also interesting to inspect whether KK graviton exchange has any influence on 
the particle densities at ground level. A summary of the ground lateral distributions 
of proton showers at vertical incidence = 3 x 10^° eV) is reported in Tabic 2. Here, 
the ground array was located at 875 g/ cm^ and the magnetic field was set to reproduce 
that prevailing upon the Auger experiment. The ratio between the mean density of 
charged particles p^j^/ p^h^"^ is a monotonically decreasing function of the distance to 
the shower core R. Nevertheless, one should note that within the error hmits the ratio, 
py pf^^^ = A*^) 7, all charged particles), is always consistent with 1. Thus, we 
deduce that the possible signatures of KK emission are entirely hidden when the shower 
front reaches the ground. Furthermore, the competition of decay and interaction of 
the first generations of mesons propagating in a medium with varying density profile, 
is of particular relevance in the indirect analysis of data collected by ground arrays. 
As can be seen in Fig. 4, the ratio would generally depend on the primary 

zenith angle. 

Putting all this together, KK-graviton exchange offers a viable mechanism to reduce 
by around 6% the mean free path of ultra high energy (E > 5 x 10^^ eV) hadrons in 
the atmosphere. 

4 Conclusion 

Theories with large compact dimensions and TeV-scale quantum gravity represent 
a radical departure from previous fundamental particle physics. If these scenarii have 
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i? = 50 m 


= 500 m 


Pch 


269.05 ± 19.9 X 10^ 


14.41 ± 1.51 X 10^ 




138.02 ± 5.37 X 10^ 


19.43 ±0.86 X 10^ 




267.56 ± 19.9 x 10^ 


12.41 ± 1.50 X 10^ 


P7 


148.85 ± 3.90 x 10^ 


25.02 ± 1.04 X 10^ 



(4±n) 


i? = 50 m 


i? = 500 m 


Pch 


245.98 ± 7.97 x 10^ 


15.09 ± 1.58 X 10^ 


Pn± 


136.85 ±6.26 x 10^ 


19.43 ±0.90 X 10^ 




244.49 ± 7.93 x 10^ 


13.11 ± 1.56 X 10^ 


Pi 


150.96 ± 4.46 x 10^ 


27.10 ± 1.94 X 10^ 
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Figure 4: Atmospheric cascade developments of proton showers for extreme primary 
zenith angles (0° and 50°) and E' = 3 x 10^'' eV. The error bars indicate RMS fluctua- 
tions of the means. 
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some truth, the scattering phenomenology above coUider energies would be quite dis- 
tinct from SM expectations. In particular, the exchange of KK towers of gravitons 
leads to a modification of SM hadronic cross sections at s^^^ > 100 TeV. Extremely 
high energy cosmic rays that impinge on stationary nucleons at the top of the atmo- 
sphere start chain reactions where the cm. energy can be as high as 500 TeV. It is 
therefore instructive to explore KK exchange sensitivity within the entire average pro- 
file of the air shower. In this paper we have contributed a few results to this question. 
We have shown that the exchange of KK gravitons could affect the rate of development 
of atmospheric cascades initiated by protons. For primary energies above 3 x 10^° eV, 
the effects are statistically significant and can thus be observed by fluorescence de- 



tectors ||3T|. We have also proved that the footprints left by the (4 + n)-dimensional 
gravitons become washed out as the shower front gets closer to the ground and, in gen- 
eral, cannot be traced back with surface array data. The details of our analysis should 
be treated with some caution since they may be sensitive to the hadronic interaction 
model used. The overall conclusion, however, should remain the same. 
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